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We show that galaxy redshift surveys sensitively probe the neutrino mass, with eV mass neutrinos 
suppressing power by a factor of two. The Sloan Digital Sky Survey can potentially detect N nearly 
degenerate massive neutrino species with mass m„ > 0.65(f2 m /i 2 /0.1./V) ' 8 eV at better than 2a once 
microwave background experiments measure two other cosmological parameters. Significant overlap 
exists between this region and that implied by the LSND experiment, and even m v ~ 0.01 — 0.1 eV, 
as implied by the atmospheric anomaly, can affect cosmological measurements. 
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Current neutrino experiments reveal anomalies resolv- 
able by nonzero neutrino masses and flavor oscillations. 
The LSND direct detection experiment suggests v^ to v e 
oscillations with <5m 2 ^ 0.2 eV 2 Q. The deficit of \x neu- 
trinos in atmospheric showers indicates mixing between 
Vp and another species with 5m 2 ^ ~ 10~ 3 — 10~ 2 eV 2 
Mi. Finally, the solar neutrino deficit requires 5m 2 ei ~ 
10 _5 eV 2 (see e.g. H for recent assessments). These re- 
sults are consistent with one to three weakly interacting 
neutrinos in the eV mass range 0). 

Cosmological measurements provide an independent, 
albeit indirect |5| , means of determining neutrino masses 
in the above-mentioned range. Massive neutrinos would 
produce a strong suppression in the clustering of galaxies, 
with even a 10% neutrino contribution making a 100% 
difference in the power p|. Detecting this suppression 
would measure the absolute mass of the neutrinos, in 
contrast to the mass splittings measured by the oscilla- 
tion effects described above. 

While the general effect is well known, most work to 
date has focused on a combined neutrino mass around 
5 eV, as this is the minimum needed to affect cosmol- 
ogy at the current observational sensitivities MM . There 
are three reasons why this situation is likely to change 
soon. First, evidence continues to mount that we live 
in a low-density universe (e.g. 0). Since the cosmolog- 
ical effects depend on the density fraction supplied by 
neutrinos, our sensitivity to the neutrino mass increases 
roughly in inverse proportion to the density parameter. 
Second, the cosmic microwave background (CMB) ex- 
periments currently under development should establish 
a cosmological framework (e.g. jHJ) that is as secure as 
the standard model of particle physics. The parameters 
left unspecified by the model may then be measured with 
confidence. Finally, upcoming high precision galaxy sur- 
veys such as the Sloan Digital Sky Survey (SDSS) |ll]] 
should be able to measure the total power on the rele- 
vant scales to ~1% accuracy. The combination of these 
developments implies that galaxy surveys will soon pro- 
vide either an interesting constraint on or a detection of 
the mass of the neutrinos. 

Although the effects of massive neutrinos are very 
large, variations in other cosmological parameters may 



mimic the signal. Therefore, to qualify as a true detec- 
tion, all other aspects of the cosmology that similarly 
affect the power spectrum must be previously or simul- 
taneously determined. 

In this Letter, we evaluate the ability of galaxy surveys 
to distinguish between these possibilities and thereby 
measure the mass of the neutrinos. We establish the 
physical basis of the measurement, evaluate the uncer- 
tainties caused by our ignorance of other aspects of cos- 
mology, and present the 2a detection threshold in mass 
for SDSS. These results depend on two assumptions: 
that CMB observations will confirm that structure forms 
through the gravitational instability of cold dark matter, 
and that the galaxy bias is linear, i.e., the galaxy power 
spectrum is proportional to the underlying mass power 
spectrum. The second assumption is relaxed in the con- 
cluding remarks. 

Neutrino Signature. -- In a universe with the stan- 
dard thermal history |13] , the temperature of the back- 
ground neutrinos is (4/11) 1 / 3 that of the CMB. This im- 
plies Oj,/i 2 sw J ?Vm„/94eV for N massive neutrino species 
of nearly identical mass to„. Here and below, Qi is the 
fraction of the critical density contributed by the zth mat- 
ter species (v = neutrinos, b — baryons, m = all matter 
species) and H = WOh km s _1 Mpc -1 . We assume this 
thermal history and a power-law spectrum of initial adi- 
abatic density fluctuations with P(k) oc k n throughout. 

This initial power spectrum is processed by the grav- 
itational instability of the fluctuations. The large mo- 
mentum of cosmological eV mass neutrinos prevents 
them from clustering with the cold components on scales 
smaller than the neutrinos can move in a Hubble time. 
The growth of the fluctuations is therefore suppressed on 
all scales below the horizon when the neutrinos become 
nonrelativistic H 
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The small-scale suppression is given by 
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Galaxy surveys such as the SDSS Bright Red Galaxy 
(BRG) survey (assumed to be volume-limited to 
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FIG. 1. Effect of a 1 eV neutrino on the BRG power spec- 
trum compared with expected precision of the SDSS BRG 
survey (la error boxes, assuming erg = 2 for the BRGs). Up- 
per curves: an Q, m — 1.0, h — 0.5, Qi,h 2 = 0.0125, n = 1 
model with and without a 1 eV neutrino mass. Lower curves: 
the same but for an fl m — 0.2, h — 0.65 model. 



lft -1 Gpc uW) should measure the power between 0.1 — 
0.2ft, Mpc~ to ~1%. A more detailed analysis shows that 
only masses below 



0.02(O m ft 2 /0.1iV) e V 
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make less than a 2ct change in the power spectrum mea- 
sured by the BRG survey. 

As an example, we plot in Fig. R] the power spectrum 
with and without a single 1 eV massive neutrino species 
for an Q m = 0.2, ft = 0.65 model (lower curves) and an 
ft m — 1.0, ft — 0.5 model (upper curves). The expected 
lcr error boxes from the BRG survey shows that the two 
models are clearly distinguishable. For comparison, the 
difference between these models in the CMB power spec- 
trum at degree angular scales is roughly 3% (1%) and 
never exceeds 5% (4%) for multipoles t < 2000 for the 
open variant of the low (high) Ooft 2 cases (c.f. |l3|). 

Parameter Degeneracies. — Although the suppression 
of power caused by massive neutrinos is large, we must 
consider whether other cosmological effects can mimic 
this signal. The suppression begins at k nr [Eq. (jl])] 
and approaches the constant factor of Eq. (0) at smaller 
scales. Many cosmological effects can produce the gross 
effect of a change in the ratio of large to small scale power; 
we must rely on the detailed differences between these 
mechanisms in order to distinguish one from another. 

We consider variants of the adiabatic cold dark matter 
model. The power spectrum is then a function of the nor- 
malization A, tilt n, ft, £l m h 2 , f^ft 2 and f^ft 2 . The spa- 
tial curvature, cosmological constant, and the linear bias 
parameter are implicitly included through the normal- 
ization H. We estimate the accuracy with which these 
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FIG. 2. Standard deviation a(m u ) as a function of the up- 
per cutoff fc max for several different choices of prior cosmolog- 
ical constraints. Models are the same as in Fig. Py and have 
m„ = leV. (a) High fl m h 2 : no priors, solid line; single prior 
of a(Q m h 2 ) — 0.04, dashed line; full CMB prior (see text), 
long-dashed line, (b) Low Q, m h 2 : as (a), save that the single 
prior is a(n) — 0.06, dashed line. 



parameters can be jointly measured from the SDSS BRG 
survey using the technique described in J14J. Here the 
6x6 covariance matrix of the 6 parameter estimates is 
approximated by the inverse of the so-called Fisher infor- 
mation matrix F. Its elements F^ are obtained by inte- 
grating {pi In P) (dj In P) over the range k m i n < k < /c max 
discussed below, weighted by a function that incorporates 
the relevant aspects of the survey geometry and sampling 
density. Here the derivatives are with respect to the i th 
and j th parameters, evaluated at a fiducial model. 

If a small variation in one parameter can be mimicked 
by joint variations in other parameters, then one of these 
functions di In P can be approximated by a linear com- 
bination of the others. This situation is referred to as a 
parameter degeneracy, since it makes F nearly singular 
and leads to extremely poor determinations (large vari- 
ance F^ ) for the parameters involved. 

Clearly, the ability to estimate parameters comes only 
from scales on which one has both precise measurements 
and reliable theoretical predictions. On large scales, lin- 
ear perturbation theory is accurate, but the survey vol- 
ume (about l/i~ 3 Gpc 3 for the BRG survey) is limited; 
hence k m i n ~ 0.005ft Mpc -1 . On small scales, linear 
theory fails near k « 0.2ft Mpc -1 . While we expect 
that detailed data analysis will push fc max to slightly 
smaller scales by including mild corrections to linear the- 
ory, we simply use the linear power spectrum for this 
work and adopt fc max ~ 0.2ft Mpc -1 . We shall see that if 
&min ^ k nT ;$ fc max , then the unique signature of massive 
neutrinos can be identified and m^ measured. 

The solid lines in Fig. show the standard devia- 



tion of a measurement of m v (or, equivalently, fi„ft 2 ) 
as a function of fc max if all relevant cosmological pa- 
rameters are determined simultaneously from the SDSS 
BRG data set. Consider first the low fl m h 2 case (bot- 
tom panel). The standard deviation drops rapidly near 
fcmax = 0.05ft Mpc -1 , well below the scale at which the 
neutrinos begin to affect the power spectrum (c.f. Fig. PA). 
If we use only information from k < 0.05ft. Mpc~ , we find 
that the neutrino signal can be accurately duplicated by 
variations in other parameters. For example, a change in 
normalization and tilt would be indistinguishable within 
the BRG survey error bars. When considering smaller 
scales, more subtle combinations still exist; these near 
degeneracies reduce the parameter sensitivity more than 
100-fold. A similar situation occurs in the high £l m h 2 
case (top panel) but at a smaller scale [Eq. (|l|)]. 

If we possess external information on the other cosmo- 
logical parameters, the situation improves dramatically 
because parameters may no longer be shifted arbitrarily 
so as to mimic the neutrino signal. Indeed, upcoming 
CMB anisotropy experiments should yield precise mea- 
surements of cosmological parameters critical to this sit- 
uation. We therefore show the effect (long-dashed curve) 
of placing CMB constraints on the cosmological param- 
eters: o-(lnA) = 0.40, a(n) = 0.06, a(Q, m h 2 ) = 0.04, 
a(£lbh 2 ) — O.lilb/i 2 , and cr(h) — 0.1, where a(i) denotes 
the standard deviation of i. We view these constraints as 
quite conservative, since they are weaker than those pre- 
dicted for the MAP satellite |lq,W and ignore the tight 
correlation between the marginalized error bars |17j . 

Which one prior is most important depends upon the 
fiducial model. For low £l m h 2 models with small neu- 
trino fractions, one cannot accurately probe the scales on 
which the neutrino suppression is small since fc m ; n > fc nr - 
This enables the tilt n to produce much of the desired 
effect. We show the error bars resulting from including 
only the tilt prior in Fig. Q (short-dash line). While this 
is the most important prior at /c max = 0.2ft, Mpc~ , the 
others combined have a non-negligible effect. 

For the high fl m h 2 case, one has precise measurements 
around k nl -, so that the onset of neutrino effects can be 
distinguished from tilt (c.f. Fig. [jj). However, altering 
J7 m ft 2 or ft, causes P(k) to slide horizontally (leaving the 
largest scales unchanged); as one can see in Fig. || this is 
roughly degenerate with the neutrino effect. The fi m ft 2 
prior is most important in this case; we show this situa- 
tion in Fig. @ (upper panel, short-dash line). 

We also test how a{m v ) increases as we double each 
prior in turn. The results change by more than a few 
percent only for £l m h 2 (20%) and n (10%) in the high 
f2 m ft 2 model and for n (40%) in the low fl m h 2 model. 

Results. - Given the confusion with variations in 
other cosmological parameters, what is the minimum de- 
tectable neutrino mass mdct? In Fig. H, we show the 
2a detection threshold [i.e. m„ = 2cr(m„)] assuming the 
CMB priors given above, fc max = 0.2ft Mpc -1 , and a fam- 
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FIG. 3. The 2a detection threshold for m„ from the SDSS 
BRG survey as a function of the matter density Q. m h 2 for 
the number of degenerate mass neutrinos N = 1-3. We have 
used h = 0.5, Q. b h 2 = 0.0125, n = 1, and fc max = 0.2/iMpc" 1 ; 
variations on these produce only mild shifts. 

ily of fiducial models with fii,h 2 — 0.0125, ft, = 0.5, and 
n = l. The choice of a fiducial model does not amount 
to fixing cosmological parameters; all parameters are de- 
termined by the galaxy data or by the prior constraints. 
With these choices, SDSS can detect the neutrino(s) if 



m u > rridct ~ 0.65 
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If the exponent here were unity, it would correspond to a 
fixed fractional suppression of power [Eq. (0)]. In prac- 
tice, one does slightly better at larger Q m hr/N because 
fc nr [Eq. (Q)] is larger and thus better resolved. 

This result is fairly insensitive to changes in the fiducial 
model or survey parameters. Choosing ft = 0.8 increases 
"idct by 15% at low fl m h 2 ; doubling the baryon density 
does the opposite. Neither matters at high Q, m h 2 . Alter- 
ing n or A affects the answer very little. Reverting from 
the deeper BRG survey to the main SDSS North survey 
|ll| increases mdet by less than 25%. 

As for the assumptions implicit in Fig. 0, only the prior 
constraints on the tilt n in the low floh 2 regime and HqH 2 
itself in the high regime are essential. We have taken 
conservative priors from the CMB here and save a full 
joint analysis for future work [n7|. Decreasing fc max to 
0.13ft Mpc" 1 increases mdct by ^40% at large Q m h 2 but 
makes little difference at small f2 m ft 2 (c.f. Fig. m). 

Quasi-linear evolution near fc max presents a complica- 
tion to the analysis, but so long as the power spectrum 
can be calculated as a function of cosmological parame- 
ters through simulations or analytic approximations p8| , 
this need not necessarily degrade the parameter estima- 
tion. However, to the extent that evolution washes out 
features in the power spectrum, degeneracies may appear 



that are not present in linear calculation. Hence, this is- 
sue merits further investigation, although we view our 
choice of fc max as conservative. Fortunately, as shown in 
Fig. 0, prior information from the CMB assists in making 
the results robust against changes in fc max . 

Our fundamental assumption is that the power spec- 
trum of the galaxies is proportional to that of the un- 
derlying mass, i.e., that the galaxy bias is linear. This 
assumption is well-motivated in the linear regime p9[ . 
Bias that develops a scale dependence as fluctuations be- 
come non-linear introduces only moderate uncertainties 
because it does not accurately mimic the neutrino signa- 
ture. Adopting the prescription of PCJ [their eq. (20)] 
and marginalizing over the scale-dependence parameter 
&2 degrades the result in the models of Fig. g by a factor 
of 2.5 implying a similar effect on the detection threshold 
of Eq. (||). The degradation stems not from the presence 
of non-linear bias but from our ignorance of its ampli- 
tude. The latter can be constrained by the scale de- 
pendence of redshift space distortions and the relative 
bias of different galaxy populations. Peacock |2]J deter- 
mines <r(&2) = 0.01 for the relative bias of IRAS and 
optical galaxies; bounding scale-dependent bias at this 
level would restore the limits of Eq. (Q) . 

In summary, although galaxy surveys can measure 
power to ~1%, isolating the mass of the neutrino at 2<r 
requires ^50% power variations. In principle, this means 
that a 50-fold improvement of mdct to m m i n [Eq. (p|)] 
would be available if other cosmological parameters were 
known perfectly. While this is unrealistic, a measure- 
ment with a precision of cr(ln A) ~ 0.1 and a(n) ~ 0.03 
would yield a factor of two improvement in rndct and is 
potentially within reach of planned experiments. 

For m v between m m ; n and TOdeti the effects on the 
power spectrum are significant yet cannot be robustly at- 
tributed to massive neutrinos. As this brackets the mass 
range implied by atmospheric neutrinos, the possibility 
of a light massive neutrino species must be considered 
when measuring other cosmological parameters. 

Massive neutrinos present an example where the 
galaxy power spectrum provides cosmological informa- 
tion on fundamental physics not available in CMB 
measurements, but where CMB measurements are are 
nonetheless needed for an unambiguous detection. This 
illustrates the complementary nature of galaxy surveys 
and CMB anisotropies. 
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